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Abstract 

Several theories for weakly damped free-surface flows have been formulated. In this 
paper we use the linear approximation to the Navier-Stokes equations to derive a 
new set of equations for potential flow which include dissipation due to viscosity. A 
viscous correction is added not only to the irrotational pressure (Bernoulli's equa- 
tion), but also to the kinematic boundary condition. The nonlinear Schrodinger 
(NLS) equation that one can derive from the new set of equations to describe the 
modulations of weakly nonlinear, weakly damped deep-water gravity waves turns 
out to be the classical damped version of the NLS equation that has been used by 
many authors without rigorous justification. 
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1 Introduction 



Even though the irrotational theory of free-surface flows can predict success- 
fully many observed wave phenomena, viscous effects cannot be neglected 
under certain circumstances. Indeed the question of dissipation in poten- 
tial flows of fluid with a free surface is a very important one. As stated by 
Longuet-Higgins I (119921 ). it would be convenient to have equations and bound- 
ary conditions of comparable simplicity as for undamped free-surface flows. 
The peculiarity here lies in the fact that the viscous term in the Navier-Stokes 
equations is identically equal to zero for a velocity deriving from a potential. 



Th e effect of v i scosity on free osc i llator y waves on deep water was st udied 
by iBoussinesq I ( 118951 ) and lLamb I (119321 ). among others. iBasset I (118881 ) also 
worked on viscous damping of water waves. It should be pointed out that the 
famous treatise on hydrodynamics by Lamb has six editions. The paragraphs 
on wave damping are not present in the first edition (1879) while they are 
present in the third edition (1906). The authors did not have access to the 
second edition, so it is possible that Boussinesq and Lamb published similar 
results at the same time. Lamb derived the decay rate of the linear wave am- 
plitude in two different ways: in § 348 of the sixth edition by a dissipation 
calculation (this is also what Boussinesq did) and in § 349 by a direct calcula- 
tion based on the linearized Navier-Stokes equations. Let a denote the wave 
amplitude, v the kinematic viscosity of the fluid and k the wavenumber of the 
decaying wave. Lamb showed that 



dot n , 9 
— = -2uk 2 a. 
dt 



XI) 



Equation (11.11) leads to the classical law for viscous decay of waves of amplitude 
a, namely a ~ exp(— 2uk 2 t). 

In order to include dissipation accurately into potential flow solutions, one 
must s o mehow take into ac c ount vorticity. In the papers by I Joseph and Wang 
(120041 ). IWang and Joseph! (120061 ) . vorticity was taken into account only in 
Bernoulli's equation, while only the potentia l component of the velocity was 
used in the kinematic condition. iTuck I (119741 ). as is commonly done in the ship 
research community, derived a single linearized free-surface condition with a 
dissip ative term, that c ombin es the kinematic and dynamic boundary condi- 

(1199 ll ) added the vort ical component in the ki nematic 



Ruvinsky et al. 



Longuet-Higgins I (Il992l ) simpli- 



tions. 

boundary condition but did not simplify it . 
fled the equations of iRuvinsky et al. I (|199ll ) by introducing a new free surface 
differing from the classical free surface by the integral over time of the vorti- 
cal component of the velocity. Below we proceed differently and show that the 
small vortical component of the velocity plays a role in the kinematic condition 
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as a dissipative term. The new resulting set of equations, which describes po- 
tential flow with dissipat i on, lea ds to a dispersion relation which corresponds 
exactly to that of lLamb I (119321 ) in the limit of small viscosity. 



2 Derivation of the new set of equations in the linear approxima- 
tion 



In order to clearly show how the new set of equations is derived, we first intro- 
duce the correction due to viscosity in the linearized equations for the potential 
flow of an incompressible fluid with a free surface (the linearization applies 
to the dynamic boundary condition expressed through Bernoulli's equation 
and the kinematic boundary condition on the free surface). We consider an 
"almost" potential flow, with its velocity v given by 

v = V<f) + V x A, (2.2) 



where A is a vector stream function. It is assumed that the vortical component 
of the velocity V x A is much less than the potential component of the velocity 

V0. 

It is important to emphasize that we can analyse the correction due to dissipa- 
tion even in the linear approximation to the Navier-Stokes equations. There- 
fore we consider a linear and "almost" potential flow in the two-dimensional 
(2D) case. Assuming for simplicity that the water depth is infinite, the equa- 
tions for the velocity v = (u, w) read 



du 1 dp ( d 2 u d 2 u\ 

dt pdx \dx 2 dz 2 J ' 
dw 1 dp ^ / d 2 w ^ d 2 w\ 

dt pdz \ dx 2 dz 2 J ' 

with the condition of flow incompressibility 

du ^ dw q ^ ^\ 

dx dz 



Using the Helmholtz decomposition for the velocity (12.21) yields 



d(f) dA y 

U ^ Z ' t) = Tx-^z- 
, . d(f) dA v 

w{x > Z > t) = d-z + ^x-> 



(2.5) 
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since in 2D there is a single component to the vector stream function, which 
we denote by A y . Equations (12. 3p combined with equation (12.41) can be written 
as follows: 



dAy 

dt 


\ dx 2 


3 2 V 

dz 2 j 


)■ 


d<j> 


p(x,z,t) 


-gz^ 


Po 


~dt 


p 


P 



(2.6) 



To determine the 'normal modes' which are periodic in respect of x with a 
prescribed wavelength 2ir /k, we assume a time-factor e~ luJt and a space-factor 
e lkx . The solutions for the potential <fi and the single component of the vector 
potential A y are then 



<f){x,z,t)=<f) e iikx -" t) e lklz , 

fly{X, Z, L) — /±y„ti a 



where 



to 



z 1 i ■ 

m = k —i — . 



(2.7) 



(2.8) 



Equation (12.81) can be rewritten as 



m 



V2 



\ 



'k* + -2 + k 2 - % 



\ 



l^+^-k 2 \. 



(2.9) 



Note that in the inviscid limit, both the vector stream function A and the 
viscosity v are equal to zero. 

Let us now write down the boundary conditions along the free surface. The 
linearized kinematic boundary condition reads 



drj 
dt 



w(x, 0, t), 



[2.10) 



since the velocity is evaluated at z = in the linear approximation. It yields 
rj(x,t) = - (i\k\<Po - kA yQ ) e^-^. (2.11) 



The linearized dynamic boundary conditions read 
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dw 

p-2pv—=p at z = r)(x,t), 

( du dw\ n . . 

"(& + feJ =0 atz = a (212) 

The boundary conditions (12.101) and (12. 12[) provide two pieces of information: 
(1) the relationship between the potential <fto arid the vector stream function 



A 

^2/0 ' 



Ay ° ~ ^tf 00 - " 2 w 1+2^' (2 - 13) 



(2) the dispersion relation uj{k) 

In the limit of small viscosity vk 2 -C lu, the right-hand side of (12. 14j) vanishes 
and the following approximation for u can be derived: 



io = ±Jg\k\ -2ivk 2 . (2.15) 



Equation (12.151) leads to the classical law (II. lft for viscous decay of free waves. 
Under the same limit [yk 2 <C ou), equation (12. 13[) can be written as follows: 

*4>=2v-?*-. (2.16) 
dxdt dx 2 dz 



This equation was derived by iRuvinsky et al. I (1199ll ) . 



So far, we have only described the viscous solution to the linear Navier-Stokes 
equations with small viscosity and as said in the introduction it is a well- 
known solution. But at this point a natural question arises: Can we describe 
this flow by using only the potential part of the velocity? 

First, let us look at the kinematic condition (1 2 . 1 U j) and explicitly separate the 
potential and vortical components of the velocities in it: 

dt~dz + dx atZ "°- (2 - 17) 



The vortical part of the velocity dA y /dx can be written in different ways, 
but here we want to express it in terms of rj. (Note that we consider linear 
equations and can therefore freely express one function through another.) To 
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do this, one can use the solution for n (12.111) with the relation (12.131) between 
the constants 6q an d A yQ . Thus, the following relation can be easily derived: 

0A » - 2i \k\ 3 ( -) = -2k 2 un(x,t) = 2vpL. (2.18) 



8x 1 1 \luJ l + 2ik 2 ^ IK ' ' 8x< 



These equations are one of the main results of our paper. The kinematic 
boundary condition now can be written without the vortical component of 
the velocity: 



Let us now turn to the equation for the potential compon ent of the velocity. 
The d is sipative correction to thi s equation was discussed in I Joseph and Wang 
( 120041 ). IWang and Joseph! (120061 ) . It is simply Bernoulli's equation, except that 
the pressure must be replaced by 



P 



8w 



Po + 2pz/ 



'8 2 



+ 



d 2 A 



dz 2 dxdz J 



(2.20) 



Note that the vortical component in the pressure (12.201) is of order v 2 and can 
be neglected. So, we can write down Bernoulli's equation with dissipation 

86 8 2 6 

-5£ + M = - 2v d$ atz = °- (2 ' 21) 



It is easy to check that the two boundary conditions (I2.19P and (I2.2ip lead to 
the same dispersion relation as (12.151) . 



3 Fully nonlinear equations with dissipation 



In the previous section, we added viscous terms in the linear equations of the 
2D potential flow of a fluid with a free surface. It is clear that these terms can 
be added in the nonlinear equations for potential flow. Moreover the analysis 
can be extended to three-dimensional flows: 



86 1 ,„ , ,9 8 2 6 

8ri „ „ , 86 

+ Vn ■ V6 = -- + 2vAn at z = 0. 
at ' 8z 



(3.22) 
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It is well-known that the modulations of weakly nonlinear (undamped) gravity 
waves in deep water, with basic wave number fc and frequency uo(ko), can 
be described by the non-dissipative nonlinear Schrodinger (NLS) equation for 
the envelope A of a Stokes wavetrain (written here in the 2D case) 

^-m^-r^ A=0 - (3 - 23) 



Recall that A is the envelope of the normal canonical variable a. The Fourier 
harmonics a satisfy the following relation: 




(3.24) 



where r/k is the Fourier transform of the free-surface elevation and ipf. the 
Fourier transform of t he velocity potential evaluated on the free surface, see 
Zakharov et al. I (119921 ). What happens if one tries to derive a similar equation 
from the boundary conditions (12. 19j) and (12.211) ? Dissipation then appears 
naturally in the NLS equation in the following way: 



.dA lu d 2 A 1 
'~dt ~ 8k^~dx 2 ~ ~ 2 



u kl\A\ 2 A 



(3.25) 



It is important to recognize that equation (13.251) has been used by many au- 
thors but without convincing physical explanation. Authors either simply refer 
to previous papers or argue that it is the simplest way to include dissipation! 
Here we provide a clear derivation of the damped NLS equation. It is also 
important to point out that if vorticity effects are not included in the kine- 
matic boundary condition, then the resulting NLS equation looks much more 
complicated. 



4 Discussion 



Ruvinsky et al. I (119911 ) derived a system of three equations in the limit of 
small viscosity: one for the velocity potential, one for the surface elevation 



and one for the vortical c omponent of the velo city. lLonguet-Higgins I (Il992n 



simplified these equations. iDebnath I (11994 ) and ISpivak et al. I (120021 ) applied 
these equations to particular problems. 



Here we have shown that the vortical component of velocity can be excluded 
from the equations, thus leaving only two equations for a quasi-potential flow. 



7 



A small amount of vorticity has been incorporated in the kinematic boundary 
condition. 



It should be mentioned that there are several publications where the same 
dissipation is used in both Bernoulli's equation and in the kinematic bound- 
ary condition. Such a choice is justified when one is interested in the nu- 
me rical integrat i on of p otential flow equa t ions. A mong them one c a n refe r 



to 



Baker et al. I (119891) . iDvachenko et al. I (120031 ). iDvachenko et al. I (120041 ) 



Zakharov et al. I (120051 ). This type of dissipation was used not only for the 
si mulation of gravity waves but also for the study of capillary wave turbulence 
by iPushkarev and Zakharov I (119961 ). The addition of dissipation is sometimes 
used to satisfy the radiation condition. Another question of interest is the 
value that one is supposed to use for the kinematic viscosity v. The molecular 
value v = 10~ 6 m 2 /s is far too small to be significant for normal waves. In fact 
it seems that values which are of the order of 10~ 4 or even 10 -3 are commonly 
used. Therefore it is the eddy viscosity rather than the molecular viscosity 
that ought to be used for practical applications. 

Finally, the set of dissipative equations we derived can be extended to inter- 
facial waves and more generally to multi-layer configurations. 
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